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We study anomalous heat conduction and anomalous diffusion in low dimensional systems ranging 
from nonlinear lattices, single walled carbon nanotubes, to billiard gas channels. We find that in all 
discussed systems, the anomalous heat conductivity can be connected with the anomalous diffusion, 
namely, if energy diffusion is <r 2 (t) = (Aa; 2 ) = 2Dt a (0 < a < 2), then the thermal conductivity can 
be expressed in terms of the system size L as n = cLP with (3 = 2 — 2/a. This result predicts that 
a normal diffusion (a = 1) implies a normal heat conduction obeying the Fourier law (/3 = 0), a 
superdiffusion (a > 1) implies an anomalous heat conduction with a divergent thermal conductivity 
(P > 0), and more interestingly, a subdiffusion (a < 1) implies an anomalous heat conduction with 
a convergent thermal conductivity (/3 < 0), consequently, the system is a thermal insulator in the 
thermodynamic limit. Existing numerical data support our theoretical prediction. 
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"-ft seems there is no problem in modern physics for 
which there are on record as many false starts, and as 
many theories which overlook some essential feature, as 
in the problem of the thermal conductivity of nonconduct- 
ing crystals"^. This statement is still true even if 
half century has past. In fact, many questions and 
many puzzles remain open, for example, what is 
a necessary and sufficient condition for heat con- 
duction in low dimensional system to obey the 
Fourier law is still not clear. Nevertheless, some 
encouraging progress have been achieved through 
computational simulations and analytical studies 
in recent years. We now understand that in a mo- 
mentum conserved one dimensional system, such 
as the Fermi-Pasta-Ulam model, the thermal con- 
ductivity diverges with the system size, L, as -L 2 / 5 , 
while it is L 1 / 3 if the coupling with the trans- 
verse mode is present. A connection between the 
anomalous heat conduction and anomalous diffu- 
sion has been found. Some interesting models for 
control heat flow such as thermal rectifier, ther- 
mal diode, and thermal transistor have been pro- 
posed and studied. 



I. INTRODUCTION 

Heat conduction in low dimensional systems has been 
attracting increasing attentions in the last decade (see 
Review and the references therein.) The problem is of 
great interests not only from fundamental point of view 
but also from application point of view. On the one 
hand, it is still an open question that under what dy- 
namical condition(s) will the heat conduction obey the 
Fourier law. It was suggested that the exponential in- 
stability might be the necessary condition 3] , and this 
conjecture has been further confirmed by the study of 
heat conduction in a quasi one dimensional Lorenz gas 



channel However, the studies on disordered Ehren- 
fest gas channel 0, triangle gas channel 6], polygonal gas 
channelQ, and an alternate mass hard-core potential 
chainQ show that a system with zero Lyapunov expo- 
nent can also display a normal diffusion and the heat 
conduction obeys the Fourier law. 

Unlike the billiard gas channel, we still don't have a 
clear answer to this question in lattice models. It has long 
been believed that the nonintegrability might be a suffi- 
cient condition for the Fourier law. However, the st udy 
on the Fermi-Pasta-Ulam (FPU) model0,[I3OIIll3 
shows that the heat conductivity diverges with system 
size L as (fl > 0), which means that the heat con- 
duction violates the Fourier law. It is clear now that 
the anomalous thermal conductivity is due to the soli- 
tary wave propagation along the chain [Tl|. It has been 
shown theoretically from the mode-coupling theory [l^ 
and the Peierls equation |15| that the divergent expo- 
nent /3 should be 2/5, which is supported by the numeri- 
cal results from different groups [uj 0. For a generic 
one dimensional lattice, we now understand from the 
conclusion given by Prosen and Campbell^fJ that the 
momentum conservation will lead to a divergent thermal 
conductivity. 

The 2/5 power law divergence for one dimensional non- 
linear lattice looks like quite general. However, an anal- 
ysis from the renormalization group argue that the di- 
vergent exponent should be 1 /3 for a momentum con- 
served one dimensional system [l7j. This looks like a puz- 
zle at first glance. However, if we look carefully at the 
analysis jlj}, we find that, a random noise added to each 
particle of the lattice is not possible in a one dimensional 
system. Recently, Wang and Li[l|| clarified this puzzle 
via mode-coupling theory by studying a quasi Id systems 
by allowing the lattice have both longitudinal and trans- 
verse motions. According to Wang and Li's theory, the 
1/3 power law comes from the mode coupling between 
the longitudinal modes and transverse modes. When the 
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coupling is absent, like in the ID FPU model, the 2/5 is 
recovered. 

If the momentum conservation is broken by introduc- 
ing an on-site potential like the Frenkel-Kontorova (FK) 
modeiyjj and the 4 model[T2|,|2(|, the heat conduction 
then obeys the Fourier law. 

On the other hand, with the fast development of nan- 
otechnology, many low dimensional systems such as the 
nanotubes, nanowires and nanobelt etc can be fabricated 
in the laboratory. It is necessary to know the properties 
such as the mechanical properties, electronic properties, 
and thermal properties before they can be put into any 
practical use. However, compared with other two prop- 
erties, much less is known about the thermal ones|2ll|. 
In fact, the very basic question such as whether the heat 
conduction in such systems can obey the Fourier law is 
still not clear, even th oug h an increasing studies have 
been devoted toUHH Hi El- 

Moreover, the understanding of the underlying heat 
conduction mechanism will shed lights on the applica- 
tion for designing novel thermal devices. Indeed, some 
promising works have been done in this direction. The 
first model for thermal rectifier was proposed by Terrano, 
Peyrard, and Casati26j. This model allows heat flows 
anisotropically along the chain. Most recently, based 
on two coupled FK chains, a thermal diode model is 
proposed|23. This model acts like a conductor when 
a temperature gradient is applied, while it acts like an 
insulator when the gradient is reversed. In addition, in 
a certain parameter range, a negative differential ther- 
mal resistance is observed. Based on this property, a 
thermal transistor|28| model is introduced to control the 
heat flow by adjusting either temperature at or current 
through the third terminal, just like we do for charge flow 
in a MOSFET (Metal-Oxide-Semiconductor Field Effect 
Transistor) . 

In this paper, we give a detailed discussion on the 
anomalous heat conduction and anomalous diffusion in 
different systems: nonlinear lattices, single walled nan- 
otubes, and billiard gas channels. The main analytical 
result, Eq. ©, for the relation between anomalous heat 
conduction and anomalous diffusion has been reported 
in Ref29]. The paper is organized as the follows. In 
Sec II, we give a general relation between anomalous dif- 
fusion and anomalous heat conduction. In Sec III, we 
demonstrate the heat conduction and diffusion in several 
nonlinear lattices such as the FPU model and the FK 
model etc.. In Sec. IV, we discuss the heat conduction 
in quasi one dimensional lattice with transverse motion 
and the anomalous diffusion in a single walled nanotubes 
(SWNT). In Sec V, we give a detailed discussion on the 
heat conduction and diffusion in several billiard gas chan- 
nels. In the last part, Sec VI, we give conclusions. 



II. ANOMALOUS DIFFUSION AND 
ANOMALOUS HEAT CONDUCTION 

The conclusion made by Prosen and Campbell 16] 
states that in a momentum conserved one dimensional 
system, the thermal conductivity diverges with the sys- 
tem size L. However, the theory does not give any indi- 
cation that how the thermal conductivity diverges with 
system size L. 

Here we would like to find a microscopic origin of the 
anomalous heat conduction observed in many ID models. 
We shall not restrict to any specific model. This should 
give us a more general way to understand the heat con- 
duction in ID systems. 

To establish a connection between the microscopic pro- 
cess and the macroscopic heat conduction, let's consider 
a ID model of length L whose two ends are put into 
contact with thermal baths of temperature Tj, and Tji 
for the left end and the right end, respectively. Suppose 
the energy is transported by energy carriers, phonons in 
lattices and particles in billiard gas channels, from left 
heat bath to the right heat bath and vice versa. The 
mean square of displacement of the carrier, with velocity 
v, inside the system can be described by, 

(Ax 2 ) = 2Dv a t a . (1) 

The average time for the particle travels from the left end 
to the right end is 

(t LR ) « (t RL ) ~ i 2/Q - (2) 

We should stress that, the time given above is the av- 
eraged first passage time [23, |3(|. It is taken over the 
distribution of the first passage time, or a truncated dis- 
tribution of the first passage time if the average is not 
convergent like in the case of subdiffusion[30( (see also 
Ref. [3l|]) for more regorous discussion.). Moreover, time 
(J2J) is not the recurrence time because when the parti- 
cle arrives at the opposite end its motion will be "inter- 
rupted" due to of the interaction with the heat bath. It 
should be pointed out that in the cases we studied in this 
paper the exponent a in Eq. and are same. 

The heat current, i.e. the energy change between the 
two heat baths in unit time, for a single particle moves 
from left to right and comes back is: 

. = E L - E R T L -T R = _AT_ 

JS {tLR) + (tRL) " (tLR) (tLR)' 1 ' 

For L particles in the system, the heat current is, 

j = Lj s oc -ATL 1 - 2 ^, (4) 
then the heat conductivity is given by 

k = -j/VT cx L 13 (5) 

with, 

/? = 2-2/<x (6) 
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This relationship can be derived in a more rigorous 
way by using the fractional derivative |29j. It connects 
heat conduction and diffusion quantitatively. The main 
conclusion is that an anomalous diffusion indicates an 
anomalous heat conduction with a divergent (convergent) 
thermal conductivity. More precisely, our result tells us 
that: a ballistic motion means thermal conductivity pro- 
portional to the system size L, a normal diffusion means 
a normal heat conduction obeying the Fourier law, a su- 
pcrdiffusion means a divergent thermal conductivity, a 
subdiffusion means a zero thermal conductivity in the 
thermodynamic limit. 



As the derivations of above equations are based on cer- 
tain conditions, the relation Eq. JSJ is not expected to 
apply to all anomalous transports. We believe that in 
other cases, where a in Eq. (J2J) is not the same as that 
one in Eq. ifTJl. a. — 8 formula may be different from Eq. 
©. In this sense, a in Eq.Q is the same a in Eq. Q can 
be regarded as our pre-assumption, and our motivation is 
to investigate, under this assumption, what a connection 
between a and (3 can be set up. The assumption that 
the a in both Eq. JQ) and Eq. (|2J) are the same one is 
based on the fractional differential equation's description 
of the diffusive process [33, which is found to be useful in 
describing the anomalous transport in both chaotic and 
quasi-chaotic systems. However, to what extend, and 
how the fractional differential equations can be applied 
is an interested topic and is still far from complete un- 
derstanding. Our studies in Ref.pij] and in this paper 
showed that at least in the examples we have studied (in 
certain parameter regimes), the prediction based on the 
fractional differential equation is in good agreement with 
the numerical simulations. 



III. ANOMALOUS HEAT CONDUCTION AND 
DIFFUSION IN NONLINEAR LATTICE MODELS 

To study diffusion in a lattice model, we first thermal- 
ize the system to an equilibrium state with temperature 
T (corresponds to energy E ), then the middle particle is 
given a much higher temperature Tg (higher kinetic en- 
ergy) . The evolution of the energy profile along the chain 
is then recorded afterwards. Quantitatively, the width of 
the pulse can be measured by its second moment 

2 _ J(E{x,t) - E )(x - x ) 2 dx 
a W f(E(x,t)-E )dx ' [ > 

where E(x, t) is the energy distribution at time t. xq 
is the position of energy pulse at t — 0. Generally, the 
energy profile spreads as 

<7 2 {t) = 2Dt a , (0<a<2). (8) 

To suppress statistical fluctuations, an ensemble aver- 
age over many realizations is performed when obtaining 
a 2 (t). 

Generally, one cannot define the same quantity for en- 
ergy diffusion in lattice models as Eq.QJ like we do in 
billiard gas channels. In billiard gas channel, as we will 
show later on that (Ax 2 ) is calculated over the particle 
concentration while a 2 in Eq. Q) is calculated over the 
energy density distribution. In fact, in billiard gas chan- 
nel, as the particle is the energy carrier, it can be shown 
that both (Ax 2 ) and cr 2 (t) are proportional to t a . How- 
ever, for lattice model, the physical picture for the (Ax 2 ) 
is not clear, thus we use Eq. Q) to describe the diffusion. 

A. The Fermi-Pasta-Ulam (3 model 



It is worth mentioning that, in past decades, a large 
amount of efforts have been devoted to the study of 
anomalous diffusion and its relationship with anoma- 
lous electric conductivity (see Ref. [33| and the refer- 
ence therein). So far, our workj^j has been the first one 
to connect the anomalous diffusion with anomalous heat 
conduction. More importantly, recent years have wit- 
nessed renaissance of the research on anomalous diffusion 
and its connection with underlying dynamics. For more 
information, please refer to a recent extensive review [3l| 
and the references therein. 



In the following section, we investigate the diffusion 
in different nonlinear lattice models, in particular, the 
FPU model and the FK model, two representatives for 
nonlinear lattices without and with on-site potential, to 
see how the anomalous/normal diffusion is connected to 
the anomalous/normal conductivity. We also check this 
relation in single walled nanotubes in Section IV and in 
the billiard gas channels in Section V. 



As is well known that the FPU model plays an im- 
portant role in the development of nonlinear science and 
the development of molecular dynamics (see the other 
articles in this Focus issue). We should say that, the 
FPU model plays also an essential role in initiating the 
study of heat conduction in low dimensional systems in 
last years. In fact, at the start stage of the study in re- 
cent years, most works have been focused on the FPU 
modelH Ea [H 13 13 III El The Hamiltonian of the 
Fermi-Pasta-Ulam fj model is 

H = y~] — - + -(x i -x i - 1 -a) 2 + — (x i -x i - 1 -a)' 1 (9) 

or 

H = E ^ + - *-i) a + t (?< " 9, - l)4 ' (10) 

i 

if the relative coordinates, <&(= Xi — id), are used. 

The heat conduction has been studied by different 
groups [9l IToL ITU fl2l fl3| numerically and found that k ~ 
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iV°- 42 ~ - 45 . Both mode-coupling theory from Lepri^jj 
and the Peierls equation study by Perezenev^ argued 
that the divergent exponent should be 2/5 which is very 
close to the numerical values. Here we give at the first 
time the diffusion behavior of the FPU model and try to 
testify the formula Eq. © in the FPU model. 

o 2 (t) versus t for different temperatures are shown in 
FigH It is shown that the diffusion takes a smooth 
crossover from ballistic one (with slope=2) at low tem- 
perature to an anomalous diffusion (slope =1.25) at 
high temperature. It is noted that at high temperature 
regimes (T > 3.0), the slope 1.25 does not depend on the 
temperature. It seems that this 1.25 is the high tempera- 
ture limit. The slope 2 is the low temperature limit which 
is easy to understand as at very low temperature, the an- 
harmonic term in the FPU model can be neglected. We 
would like to point out that in doing the best fit for the 
slopes given in Fig^ one needs to get rid of the bound- 
ary effects (same is true for Figs. 3, 5 and 9). This can be 
easily down by observing the time when the front of the 
energy wave packet approaches the boundary. This time 
is used as the cutoff time for the best fit. 

In fact, it is only a qualitative way to check when 
a wavepacket front approaches the boundary from the 
snapshot. The quantitative way is to use c 2 (i). In the 
log-log scale <J 2 (t) plot, the a 2 (t) increases with time t 
until a critical time when the wave front hit the bound- 
ary and the reflection happens. After this critical time, 
the a 2 (t) becomes saturate, some time even oscillating. 

If we insert a = 1.25 into formula ©, we obtain 
(3 = 0.4 which is exactly the same value obtained 
by mode-coupling method[Tj| and the Peierls equation 
method[l5j. and very close to the numerical value, (3 = 
0.42 ~ 0.45, from direct nonequilibrium thermal conduc- 
tivity calculation. 

Snapshots of two representative pulse propagation 
along the FPU lattice are shown in Fig. Eta) for low 
temperature and Fig. \^h) for high temperature. The 
difference of the diffusion in two different temperatures 
are clearly seen. 



B. The Frenkel-Kontorova model 

The FPU model discussed above is a representative 
lattice model without on-site potential. The Frenkel- 
Kontorova model, on the other hand, is a representative 
for a general one dimensional lattice with on-site poten- 
tial. It describes a chain of particles connected by spring 
of constant 7, and put on the substrate. The average par- 
ticle distance is a, and the periodic of on-site potential is 
b. The Hamiltonian of the standard FK model is fl9l 




FIG. 1: u 2 (t) versus time, t, for the FPU model at different 
temperatures. The temperature ranges from very low to very 
high. The chain length N = 501. 10 4 realizations are taken 
for the average. T5 = 10T. It is clearly show the transition 
from ballistic motion (slope =2) at very low temperature to 
the anomalous diffusion at large temperature (slope =1.25). 
[3* = 1. Note that the first few data for T — 0.1 deviate from 
t 2 because the large fluctuations due to thermalization. 
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FIG. 2: Snapshots of the heat pulse spreading in the FPU 
model for different temperatures, (a) T = 0.01; (b)T = 10. 
In both cases, the chain is first thermalized to a state of tem- 
perature T = 0.01 and T = 10, then a heat pulse (T s = 5T) 
is given to the middle particle. The chain length is N = 501, 
and /3* = 1. 
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H 



E 



2m 



;( x i 



V 



■cos2irxi. (12) 



Z-^ 9m 9 V 



X 



i-r 



-a) 2 - 



A 2irXi 

COS — : 



(27T) 



(11) 



Here X t is position of the i'th particle. For convenience 
of numerical calculations, the above Hamiltonian can be 



After the scaling, the period of the on-site potential 
is unity, /1 is the winding number, fj, = 1, m = 1, 
and the lattice constant is also taken to be unity in this 
paper. Of course, we can also express the Hamiltonian 
(lll|l in terms of the displacement from the equilibrium, 
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FIG. 3: o 2 (t) versus time, t, for the FK model at different 
coupling constants. The system is thermalized at first at 
T = 0.25. The chain length N = 201. 5 x I0 5 realizations 
are taken for the average. At V = 1, the energy transport 
superdiffusively, while at V = 5 the energy transports diffu- 
sively which supports the results for normal heat conduction. 



qi. Substituting Xi = ifi + qi into eq 
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In Fig. |3]we show a 2 {t) versus t in log- log scale such 
that the slope of the curve gives the value of a. The 
temperature in both cases is T = 0.25. It is clearly 
seen that the asymptotic slopes show different behav- 
iors. In the case of large coupling constant, V = 5, 
the energy transports diffusively, the heat conductivity 
is a system-size independent constant, as we observed 
before[19j. However in the case of small coupling con- 
stant, V = 1, cr 2 (i) ~ t 1 - 7 , the energy transports superdi- 
fussively. According to Eq.ljSJ, this corresponds to a di- 
vergent thermal conductivity which agrees qualitatively 
with the results from Savin and Gendelman|34j 

The snapshots also show complete different story for 
V = 1 and V = 5 cases. In the former case, the heat pulse 
splits into two pulses and propagates along two opposite 
directions. However, in the case V = 5 the pulse never 
splits, it only becomes fatter and fatter. 

As is well known that the ground state of the Frcnkcl- 
Kontoroval model can be cast into a standard map. Then 
there exists a one-to-one relationship between these two 
models (see review article[35|). For example, an incom- 
mensurate state corresponds to an invariant curve with 
an irrational winding number in the phase space and a 
commensurate state corresponds to an invariant curve 
with rational winding number. In the standard map, the 
anomalous diffusion in both x and p directions has been 
extensively studied by Zaslavsky et al 36] . It is still not 
clear and deserve further investigation how the anoma- 
lous diffusion in the standard map can be connected to 
the anomalous diffusion in the FK model observed in 
Fig.Ct). 
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FIG. 4: Snapshots of the heat pulse spreading in the FK 
model at different coupling constant, (a) V = 1; (b)V = 5. 
The other conditions are the same as that ones in Fig|3] 




FIG. 5: The<7 2 (£) versus time, t, for the 4> 4 model at different 



temperature with (3* 
asymptotically. 



fO. In all three cases, a (t) 



C. The 4 model 

Like the FK model, the <fi 4 model has an on-site po- 
tential, its (dimensionless) Hamiltonian isfT^. |2G| 



\2 i P 4 



(13) 



Due to the presence of the on-site potential that breaks 
down the momentum conservation, the heat conduc- 
tion in this model has been shown to obey the Fourier 
lawfl^. l2p| . We show here <r 2 (t) versus t for different 
temperatures in Fig [S] The slopes in all three cases are 
approximately one, thus the vibration energy transports 
diffusively, that supports the conclusion of finite thermal 
conductivity. 

For comparison, we summarize the JN (proportional 
to thermal conductivity) for different models in FigEl 
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FIG. 6: JN (proportional to thermal conductivity k) versus 
system size N for a harmonic chain (★), the FPU model (•), 
the (j) 4 model (□) and the FK model (A). The JN diverges as 
N in the harmonic model and diverges as JV ' 42 in the FPU 
model, whereas it is constant in the <f> 4 model and the FK 
model. The data for the FK model are taken from Ref . |l9|l . 
while the others from R,ef|l2l|. 

D. Coupled rotator model 

Another interesting lattice model is the so-called cou- 
pled rotator model. The Hamiltonian is: 

H = V p- + (1 - cosfo - n-i - a)), (14) 
i 

where a is the lattice constant which is taken to be 2ir. 
Obviously, due to the absence of the on-site potential, 
the momentum is conserved in this model, nevertheless, 
numerical results [33, |3£| show that the heat conduction 
obeys the Fourier law. Recently, Flach (3^ argues that 
the finite thermal conductivity might be a numerical phe- 
nomenon, the thermal conductivity of this model is diver- 
gent. He conjectures that in order to observe the diver- 
gent thermal conductivity, one needs to go to very large 
system size, so large that it may exceed today's computer 
facility. In the following, we shall clarify this puzzle. 
Generally, heat flux in a ID lattice is defined as: 

I ( N N 1 \ 

J = I 22 fi&i + 2J Ya™ 1 ^ 1 I = Jl + J2 > ( 15 ) 

where fi is the interaction force between particle i and 
i + 1. In this model, /j = sin(xi — Xi+\ — a). J\ cor- 
responds to energy flux transfered from one particle to 
the other, while J 2 the energy flux due to the particle's 



movement. Usually, in other lattice models such as the 
FPU model, the FK model and etc, (J 2 ) vanishes be- 
cause the particle vibrates only around its equilibrium 
position. However for coupled rotator model, things are 
different. At lower temperature, so low that the kinetic 
energy is not large enough to make the rotator rotating, 
(J 2 ) also vanishes. However, if the temperature is large 
enough so that the rotator can overcome the energy bar- 
rier to start rotating, then, (J2) is not zero. This part 
of energy current does not go through the system from 
high temperature end to the low temperature end. If we 
use the non-equilibrium method to calculate the thermal 
conductivity, for example, by puttiny thermal bathes on 
the two ends of the lattice, then (J2) does not contribute 
to the thermal conductivity. 



When one uses the Green-Kubo formula to calculate 
the thermal conductivity, the situation becomes differ- 
ent. In Fig. we show J-J correlation contributed by J\ 
and J 2 . (J(O)J(i)) = ((J x (0) + J 2 (0))(Ji(t) + J 2 (t))) = 
(Ji(O)Ji(t)) + (J 2 (0)J 2 (t)) + (Ji(O)Ja(f)) + (Ja(O)Ji(t)). 
The ensemble average of the third and fourth terms van- 
ish because, J\ is symmetric with space reflection, while 
J 2 anti-symmetric. Thus (J(O)J(t)) = (Ji(O)Ji(t)) + 
( J 2 (0) J 2 (i)). Clear exponential decay and power-law de- 
cay can be seen for {J x {Q)J x {t)) H3 and ( J 2 (0) J 2 (t)), 
respectively. This implies that if only J\ is taken into 
account, then one has a normal heat conduction as is ob- 
served by previous numerical simulations [33, HI]- Oth- 
erwise if both are considered, then the J-J correlation 
decays slower than t , as a result, the thermal conduc- 
tivity diverges as Flach conjectured. 
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IV. LATTICE MODEL WITH TRANSVERSE 
MOTION AND SINGLE WALLED CARBON 
NANOTUBES 

A. A polymer like quasi one dimensional lattice 

To study the cross over of the divergent exponent from 
2/5 to 1/3, Wang and LijH introduced a polymer like 
chain. The chain consists of N point particles with mass 
m on a ID lattice. The particles have both longitudinal 
and transverse motions. The lattice fixes the connec- 
tivity topology such that only the neighboring particles 
interact. The Hamiltonian is given by 

^(P> r ) = Ef| + ^£(|r w - ri |-a) 2 

i i 

+ K^2_j cos 4>i, (16) 

i 

where the position vector r = (x, y) and momentum vec- 
tor p = (p x ,Py) are two-dimensional; a is lattice con- 
stant. The minimum energy state is at (ia,0) for i = 
to N— 1. If the system is restricted to yi — 0, it is essen- 
tially a ID gas with harmonic interaction. The coupling 
K r is the spring constant; K$ signifies bending or flexi- 
bility of the chain, while <pi is the bond angle formed with 
two neighbor sites, cos = — rij_i • n^, and unit vector 
iii = Arj/|Ai-i|, Art = r l+1 - n. 

The conductivity results, jN, is shown in Fig. [S] It 
is shown that at the high temperature the thermal con- 
ductivity diverges with system size L as, L 1 / 3 , whereas at 
the low temperature regime, the thermal conductivity di- 
verges with power 2/5. A detailed mode-coupling theory 
analysis^^ tells us that at high temperature the trans- 
verse modes are strongly coupled with the longitudinal 
ones while at the low temperature regime the longitudi- 
nal modes are weakly coupled or completely uncoupled 
with the transverse mode. This explains qualitatively 
and quantitatively the crossover from 2/5 to 1/3. More 
interestingly, in an intermediate regime we also observe 
the In N divergent thermal conductivity for which a the- 
ory is still lacking. 



B. Single walled carbon nanotubes 

Carbon nanotube is one of exciting nano-scale mate- 
rials discovered in recent years. It reveals many excel- 
lent mechanical, thermal and electronic properties |40|. 
However, compared with mechanical and electronic prop- 
erties, much less is known about thermal conductiv- 
ity in nanotubes. Although some researches have been 
done [22, many important and fundamental questions re- 
main unsolved. For example, does heat conduction in 
a SWNT obey the Fourier law? If it does not, then, 
how does the thermal conductivity diverge with the tube 
length L? Is heat conduction in carbon nanotubes like 




FIG. 8: jN vs N on double- logarithmic (a) and linear-log 
plot (b). The parameters (K^,Tl,Th) of the model are, set 
B: (1, 0.2, 0.4), set E: (0.3, 0.3, 0.5), set F: (1, 5, 7), set H: 
(0, 0.3, 0.5), set J: (0.05, 0.1, 0.2). All of them have K r = 1, 
mass m = 1, lattice spacing a — 2. The straight lines on F 
and E have slope 1/3, while the slope on H and J is 2/5. This 
picture is taken from Ref fl^l . 

that one in ID lattice model or 2D lattice model or 3D 
bulk material? 

The Hamiltonian of the carbon SWNT is: 

where 14, = f c (rij)[V R (rij)+bij Va^jj )] is the Tersoff em- 
pirical bond order potential, ry is the distance between 
the i'th and the j'th atom. Vij(r,j) = Aexp(— Ar*y), 
and VAiyij) = — B exp(— urij) are the repulsive and 
attractive parts of the potential, and f c {r) = 1, for 
r < R;(l +cos ^T^ ) )/2, for R < r < S and 0,r > S. 
bij are the so-called bond parameters depending on the 
bounding environment around atoms i and j, they implic- 
itly contain many-body information. A, _B, A, R and S 
are par ameters. For detailed information please refer to 
Ref.|41|. For the comparison we also study one dimen- 
sional carbon lattice with the same interaction, namely, 
Vi,i_l = fc(xi,i-i)[V R (x iti ^i) + Vy»(x< 

Molecular dynamics simulation 25] show that at low 
temperature such as 2K, for both the SWNT and the ID 
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FIG. 9: a 2 (t) versus time, (a) ID lattice at 2K\ (b) ID lattice 
at 300K; (c)(5,5) SWNT at 2K; (d) (5,5) SWNT at 300K. 10 3 
realizations are performed to do the average. 



carbon lattice, there is no temperature gradient. It re- 
sembles the ID lattice model with harmonic interaction 
potentially • This can be understood from the Taylor 
expansion of the Tersoff potential by keeping up to the 
second order term. Because at low temperature, the vi- 
brations of atoms are very small, the potential can be 
approximated by a harmonic one, and the vibration dis- 
placement in transverse direction is much smaller than 
the one along the tube axis and can be neglected. This 
result means that the phonons transport ballistically in 
the SWNTs at low temperature. 

However at room temperature, the situation changes 
dramatically. At 300K there is still no temperature gradi- 
ent in ID carbon lattice, but in the SWNT, temperature 
gradient is set up. In ID carbon lattice with Tersoff po- 
tential, the increase of temperature does not change the 
harmonic character because the strength of Tersoff po- 
tential; while in the SWNT at room temperature, there 
are anharmonic terms due to the transverse vibrations. 

In Fig. 0we show (J 2 it) versus t in double logarithmic 
scale, so that the slope of the curve before the turning 
point (due to boundary reflection) gives the value of a. 
FigUJa) and (b) show the results for the ID lattice at 
temperature 2K and 300K, respectively. It is clearly seen 
that the slopes are 2 in both cases. This means that en- 
ergy transports ballistically. However, for the SWNT, the 
situation is different. At T=2K, energy transports ballis- 
tically (Fig. Efc)) like in the Id carbon lattice, whereas at 
T=300K, energy transports superdiffusively with a ~ 1.2 
(FigEJd)). If we use a ~ 1.2, then, according to our for- 
mula Eq. ©, we have /3 ~ 1/3 which agrees very well 
with the result, (3 » 0.32 obtained by Maruyama[23j with 
non-equilibrium molecular dynamics method. 

Snapshots of two representative pulse propagation 
along the tube/lattice are shown in Fig^| In ID carbon 
lattice, as demonstrated by Fi dlOf a). both at T=2K and 
T=300K, a peak splits into two peaks very quickly, and 
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FIG. 10: Snapshots of heat pulse spreading in a ID carbon 
lattice and a (5,5) SWNT at different temperatures, (a) ID 
lattice at 2K. (b) The SWNT at 300A". For ID lattice, at 
T=300K, the pulse spreads in the same manner and same 
speed as for T=2K, whereas for SWNT, the pulse spreads 
faster at T=300K than at T=2K. For quantitative compari- 
son, please refer to Fig 



each peak spreads ballistically to the ends of the lattice. 
However, in the SWNT, the pulse spreads out, but with 
different rate at 2K and 300K. Quantitatively, the width 
of the pulse can be measured by its second moment as is 
given by Fig El 

It is worth pointing out that, in many experiments 
or practical applications, the SWNT/lattice is put on 
a substrate, then one should add one additional poten- 
tial term in Hamiltonian (|17J) . With introduction of the 
on-site potential, heat conduction and energy diffusion 
are expected to be normal like that one in the ID FK 
model[I3 and 4 modelHHIj. 



V. BILLIARD GAS CHANNELS 

In this section, we give a brief review on recent years' 
studies on heat conduction in billiard gas channels and 
show our study on the diffusion. In fact, in addition to 
heat conduction problem^, 0, 0, 0, E|, billiard models 
have been used to study other macroscopic transport |45l 
\M . iTT l lil . l49H5fil l . For example, the Lorentz gas model 
has been used for the study of connection between the 
electrical conductivity and diffusion 01 . Moreover, in 
recent years, extensive studies on the chaotic dynamics 
and the origin of statistical law has been done by using 
different types of billiard gas model [sil Ifilj . 

In the study of heat conduction, the two ends of the 
billiard gas channel (see Fi dll|) are put into contact with 
heat baths. The heat baths are modeled by stochastic 
kernels of Gaussian type, namely, the probability dis- 
tribution of velocities for particles coming out from the 
baths is 



x \V X \ I VZ 



(18) 



for v x and v y , respectively. 

The temperature field at a stationary state is cal- 
culated by time averages by dividing the configuration 
space into a set of boxes Ci 3 . The time spent within a 
box in the jth visit is denoted by tj and the total number 
of crossings of a box Cj during the simulation is M. The 
temperature field is defined by0 



Tr 



(E)c 



(19) 



Then it is projected on x direction (the transport direc- 
tion). Here Ej(Ci) is the kinetic energy of the particle at 
its j'th crossing of the Ci box. The heat flux is calculated 
by the change of energy carried through to the left and 
right ends by the particles, 
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l M 



(20) 



where AEj = (E in — E out )j is the energy change at the 
jth collision with a heat bath, tm is the total time spent 
for M such collisions. 

In numerical simulation, we compute the flux for a 
single particle J\. The scaled heat flux is Jn(N) = 
NJi(N) 41, where N is the number of the cells. Each 
cell has length a, thus the channel has length L = Na. 

The diffusive property of the particles in the channel is 
quantified by the mean square displacement (Ax 2 (t)) — 
((x(t) — x(0)) 2 }. Here (• • • ) means the ensemble average 
over many particles. 

In the following we discuss the heat conduction and 
diffusion properties in the Lorentz gas channelH, the 
Ehrenfest gas channel 0, the triangle channel jg, the 
polygonal billiard channels, and the alternate mass 
hard-core potential modeiy|, respectively. 



(g) 



FIG. 11: The geometrical configuration of billiard gas chan- 
nels, (a) The Lorentz gas channel^- (b) The disordered 
(height) Ehrenfest channel^]- (c) The disordered (position) 
Ehrenfest channel^- (d) The triangle gas channel)^], (e) The 
triangle-square channel, (f) The polygonal gas channel^]], (g) 
The alternate mass hard-core potential model 8]. 



The dynamics in the Loretnz gas is mixing and all tra- 
jectories with nonzero projection on the x direction are 
of hyperbolic type; further it has positive Kolmogorov- 
Sinai entropy and a well defined diffusion^^ • If has been 
demonstrated Q that the heat conduction in this model 
obeys the Fourier law. 



B. The Ehrenfest gas channel 



A. The Lorentz gas channel 

In order to answer the question of whether the ex- 
ponential instability plays a crucial role in the Fourier 
law, Alonso et al|4( studied the heat conduction in the 
Lorentz gas channel (see Fig II If a) for its configuration). 
The model consists of two parallel lines of length L at dis- 
tance H, that is taken to be unity, and a series of semi- 
circles of radius R placed in a triangular lattice along 
the channel. By construction no particle can move along 
the horizontal direction without colliding with the disks. 



Instead of the semicircles of radius R, we place right 
triangles along the chain, namely, in each cell, we have 
two triangles, one on the bottom wall, the other on the 
top. The triangles are placed at the position of x = 
1,3,- • • (arbitrary unit). The height of channel is H — 
1.1, while the height of the triangle is 0.6. The unit cell 
length is 4. The channel of length N is N repetitions of 
the cell. Due to the existence of free path, this model 
demonstrate a superdiffusive motion, i.e. a 2 = (Ax 2 ) tx 
2Dt a , with a w 1.672 ±0.003 @. The heat conductivity 
diverges with the system size as k oc iV 814 . This (3 — 
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FIG. 12: Energy diffusion in the Ehrenfest gas channel, (a) 
P(x,t) versus x for a fixed t — 45, 000. Inset of (a), we show 
the scaled PDF in logarithmic scale, — log(P(x, t)/P(x, 0)). 
(b) a 2 (t) versus time t. 



0.814 is very close to that value 0.804 obtained from the 
formula eq.® by inserting a = 1.672. 

In our numerical simulations, initially the particles are 
put at x = (the middle of the channel) , with a uniform 
concentration distribution along the y direction of the 
channel. Their speeds are set to be unit with a uniformly 
distributed direction in x — y plane. Numerically, it is 
found that probability distribution function (PDF) (the 
probability of finding particle at position x at time t) can 
be described by 

P(x, t) = P(0, t) exp (- Co \x/a(t)p) , (21) 

with c = 1.58, and 7 = 0.351. In Fig. H3 
P(x,t) is drawn for t = 45,000, in the insert we draw 
— \og(P(x, t)/P(0, t)) versus x which shows a good power 
law. As we shall demonstrate later on, this formula H21|) 
seems to be quite general for all superdiffusive case in bil- 
liard channels studied in this paper. The only difference 
lies in Co, 7 and a(t). 

Channel with triangles at random positions 
In order to obtain a diffusion motion in the Ehrenfest 
gas channel, we modified above periodic structure by in- 
troducing disorder. The first way to do so is to make the 
positions of the triangles random (Fig. Ilir bl). namely, 



Xi = d * i?^, where Xi is the position away from the peri- 
odic structure. 

Channel with right triangles of random heights 

The second way is to let the heights of the triangles 
random, namely, 

hi=h + d*Ri, i = 1,2, • • • ,2N (22) 

where {i?;} are random numbers uniformly distributed in 
the interval [— 1, 1], d is the magnitude of disorder, hi < 
H, where H is the width of the channel. The geometry is 
shown in Figure lTlT c) shows the geometry of the channel. 
We take H = 1.1, h = 0.6, and d G [0,0.4]. 

Our results 5] shown that in both disordered Ehrenfest 
channels, the heat conduction obeys the Fourier law due 
to the normal diffusion induced by the disorder. 



C. The triangle gas channel 

In order to study further whether the exponential in- 
stability is a necessary condition as proposed in Ref.0, 
we consider a two dimensional billiard model which con- 
sists of two parallel lines of length L at distance H and a 
series of triangular scatterers, Fig lllf d). In this geome- 
try, no particle can move between the two reservoirs with- 
out suffering elastic collisions with the triangles. There- 
fore this model is analogous to the one studied in Lorentz 
gas channel with triangles instead of discs and the essen- 
tial difference is that in the triangular model discussed 
here the dynamical instability is linear and therefore the 
Lyapounov exponent is zero. 

The numerical results from this model depends sen- 
sitively on whether the ratio 9 /it and <p/w are rational 
or irrational numbers. In the case of irrational ratios, 
the system shows normal diffusion and the heat conduc- 
tion obeys the Fourier law, while in the case of ratio- 
nal rations, the system shows a superdiffusive behav- 
ior, (Aa; 2 ) = 2Dt 1178 as shown in Fig. 6 in Ref.@. 
The heat conductivity diverges with the system size as 
k 1=3 7V°' 25±0 01 by using larger system size. Again, 
by substituting a = 1.178 into formula ©, one gets, 
/3 = 0.30 which is close to 0.25. The deviation is due to 
the finite size effect |2flj[ . 

D. The triangle-square gas channel 

The triangle-square gas channel (see Fi dllf e') for its 
geometry) consists of isosceles rectangular triangles on 
the bottom and the rectangles on the top. The length of 
the fundamental cell is one, the height of the channel is 
H = 0.6. The height of the isosceles rectangular triangle 
is 0.4, and the width and the height of the rectangle is 
0.25 and 0.5, respectively. 

The particle's motion inside the channel is found to be 
superdiffusive. Surprisingly, we also found that the PDF 
can be described by Eq.|J2U with c = 1.95,7 = 0.376 
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FIG. 13: Energy diffusion of the triangle-square gas channel. 



and, a 2 (t) ~ f 1,34 , see Fig^l The value calculated 
from the diffusion is = 2 - 2/1.34 = 0.51. This is 
not far away from — 0.46, the direct non-equilibrium 
calculation by putting two heat baths on the two ends of 
the channel. 



E. The polygonal gas channel 

Another interesting billiard gas channel was proposed 
by Alonso et al[7|. The geometry of this model is shown 
in Fi dllf D. The height of the channel is H — 1. In 
a fundamental cell, there are two isosceles triangle with 
angle 4>2 on the bottom and a isosceles triangle with <pi 
on the top. The cell length along the horizontal direction 
is 2d, d = /i/(tan(/>i + tan 02/2). 

By fixing (f>i = (v5 — l)n/8, and adjusting the angle 
4>2 = with q = 3,4, 5, 6, 7, 8, 9, this model displays a 
very interesting transition from a superdiffusive motion 
(q = 3, 5, 6, 7), to a subdiffusive motion = 4), and nor- 
mal diffusion at q = 8 and 9. The subdiffusive motion is 
quite interesting for heat conductivity, as a < 1 in this 
case, then according to our formula ©, we immediately 
have < 0, which means that the thermal conductivity 
of such subdiffusive system is an insulator in thermody- 
namic limit. 



In the subdiffusive case, Alonso et al0 shown that a — 
0.86, and the thermal conductivity goes as k ~ L~°' 63 . 
According to our formula (©, if a = 0.86, = —0.33 
which is larger than the one obtained by Alonso et al0. 
This is not a surprise, because the channel length in their 
study of thermal conductivity is too small (L < 40). If 
the channel is longer, the value of will become much 
more closer to our theoretical estimation (0 — —0.33). 
To demonstrate this, we extend the thermal conductivity 
simulation from L e [1,40] used by Alonso et al to 
L £ [40,80], and find that = -0.48 which is more 
closer to — —0.33 than the one obtained by Alonso et 
al. If L — > oo, one can expect goes to —0.33 eventually. 



F. The alternate mass hard-core potential chain 



The billiarde gas models discussed so far do not have 
local thermal equilibrium, thus people may argue about 
the validity of temperature definition [S^]. In order to 
overcome this problem, we have recently introduced a 
gas channel that has local thermal equilibrium, and 
also exhibits normal diffusion, but with zero Lyapunov 
exponent @- 

The model is identical to the alternate mass hard-point 
gas |53l ] , namely it consists of a one-dimensional chain of 
clastically colliding particles with alternate masses m and 
M. Here, however, in order to prevent total momentum 
conservation we confine the motion of particles of mass 
M (bars) inside unit cells of size 1 = 1. Schematically 
the model is shown in Fig. Illf in which particles with 
mass m move horizontally and collide with bars of mass 
M which, besides suffering collisions with the particles, 
are elastically reflected back at the edges of their cells. In 
between collisions, particles and bars move freely. The 
total length of the system is L = Nl where N is the 
number of fundamental cells. 

Our numerical results clearly indicate that our model, 
contrary to the translationally invariant model 53] , obeys 
the Fourier law. The only difference between the two 
models is total momentum conservation. Our results sup- 
port the proof by Prosen and Campbell [Tl|. 

We shall stress an important property of diffusion 
in this model. For any m ^ M the system shows 
very nice normal Gaussian diffusion, see Fig ll4f a) for 
M = l,m = (\/5 — l)/2. However for m = M we 
found that even if a 2 (t) = 2Dt, the probability distri- 
bution function (PDF) is not a Gaussian, instead, it 
can be described by the empirical formula: AE(x, t) = 

AE{x ,t)exp (-0.826 Ix/ct^)! 1 ' 69 ) , with a 2 (t) = 2.21i. 

Except the two constants, this formula is same as that 
one for superdifussion process, Eq. I|21|) . Why these dif- 
ferent processes have similar PDF is still an open question 
which deserves further investigations. 
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FIG. 14: Energy diffusion of the alternate mass hard-core 
potential model, (a) M = 1, m = (y/E— l)/2, (b)m = M. 



VI. COMPARISON BETWEEN THEORY AND 
NUMERICAL RESULTS 

All numerical results are summarized and represented 
in Fig. 1151 where we draw (3 versus a for almost all known 
models, and compare it with Eq. JfjJ). Here we give a de- 
tail discussion by classifying the model in four categories: 
Ballistic, normal diffusion, superdiffusion, and suberdif- 
fusion. 

A ballistic motion, a = 2, leads to a divergent thermal 
conductivity k oc L. The only existing analytical result 
is heat conduction in a ID harmonic lattice. It is known 
that heat is transported by phonons in lattice model. 
Because the absence of resistance and umklapp process, 
the phonons transport ballistically in harmonic lattice 
model, thus a = 2. From our formula J^J, the thermal 
conductivity in the ID harmonic lattice diverges as 
with p = 1, this is exactly what was shown by Lebowitz 
et al|23 ("*" in Fig. ITS} . 

A normal diffusion, a = 1, means that the thermal 
conductivity is a size independent constant, (3 — 0, i.e. 
the heat conduction obeys the Fourier law. For example 
in the ID Frenkel-Kontorova model, under certain con- 




FIG. 15: The a — (3 plot. *, normal diffusion (see the text 
for more explanations). *, the ballistic transport. v/, ID 
Ehrenfest gas channel||; □, the rational triangle channel^; 
A, the polygonal billiard channel with <j>\ — (v5 — l)-7r/4) 
and 4>2 — it/3 7]; <), the triangle-square channel gas. f3 values 
are obtained from system size L G [192, 384] for all channels 
except Ehrenfest channel|f|- © the FPU lattice model at high 
temperature regime, and © the single walled nanotubes at 
room temperature. For comparison, we also draw Denisov 
et aJ's result j3 = a — 1 (dotted line) IH^| . and our results 
/3 = 2 - 2/ a (Eqinj (solid curve). 



dition (see Fig|2|), the heat energy transports diffusively, 
and the thermal conductivity is finite and independent 
of the system size L. Other ID models showing nor- 
mal diffusion and normal thermal conduction are: the 
ID Lorentz gas channel 3 , the ID disordered Ehrenfest 
gas channel 0, the ID irrational triangle channel 0, the 
alternate mass hard-core potential model||, and some 
ID polygonal billiard channels with irrationa cf>i and cer- 
tain rational triangle (f>2 0- ID disordered lattice model 
under certain boundary conditions [Hij. 'V in Fig. 1151 
represents all models with normal diffusion. 

A superdiffusion, 1 < a < 2, implies an anomalous 
heat conduction with a divergent thermal conductivity 
L 13 . The exponent < j3 < 1 differs from model to 
models. In billiard gas models, we have the Ehrenfest 
gas channel 5], the rational triangles modelQ, and the 
triangle-square channel. In lattice model, we have the 
FPU model and the nanotube model. As is shown in 
section IIIA and IVB for the FPU and the SWNT, re- 
spectively, the a and (3 values agree with Eq. © perfectly. 
(See also FiglT5*)l 

A subdiffusion. So far only one set data in this cat- 
egory is available, for which we have discussed already 
in Section VE. More models worth studying in this case 
both for billiard channels and for lattice. In fact, to find 
a lattice with subdiffusion will be an interesting problem. 
If such kind of lattices exist, they must be a thermal in- 
sulator. 
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FIG. 16: (a) The spatiotemporal PDF of the flight for the 
polygonal billiard channelfj with parameter cj>2 = 7r/3. Note 
that ifif(x,t) goes to zero for t > 1.8,which is different from 
the assumption of the power law decay in t as assumed by 
Denisov et aJ.|H^|. (b) The PDF of the first passage time 
for the same billiard with cell umber N = 4. It shows that 
4>(t) decays in a power law rather than in exponential way as 
assumed in Ref . l5al . 



A. Another scaling? 

Denisov et al(5{| has also proposed a scaling relation 
between the anomalous diffusion and anomalous heat 
conduction for the so-called dynamic channel. Accord- 
ing to their scaling, if the particle moves in a superdif- 
fusive way, then the thermal conductivity of the channel 
diverges as the system size L as, L", with 3 = a — 1. 
They state that the billiard gas channels 0,|E 00 belong 
to the class of dynamical heat channels. Here we demon- 
strate with numerical evidence that dynamical heat chan- 
nel might be useful for some unknown models, but not 
valid for the billiard gas model and lattice model. For 
comparison, we also draw 8 — a — 1 in Fig ll5l The 
deviation is clearly seen. 

In Fig HbT aL we show the spatiotemporal probability 
distribution function (PDF) of the flight (between two 
consecutive collisions with the channel boundaries) for 
the polygonal billiard channel with (f>2 = 7r/3 7J . This 
figure shows that the situation in billiard channel is more 
complicated than that described by Levy walk used in eq. 
(7) in Ref|55|. In the Levy walk, there is no correlation 
between two consecutive walks, one can thus expect a 
probability of 2~™ for a particle to move in the same 
direction (left or right) in n consecutive walks. However, 
for the polygonal billiard channel given here, we observed 
that this probability drops to zero as long as n > 6. We 
also found that the probability for a particle to 'hop' n 
cells consecutively to left (or right) deviates remarkably 
from 2~ n . This implies that the Levy walk model used 
in Ref.[5f| oversimplifies the billiard gas channel. 

As a result, there exists a fundamental difference in 
PDF of the first passage time (FPT) for Levy walk model 



and the billiard gas channels. The PDF of FPT <j>(t) for 
the same billiard model is found to decay asymptotically 
in power law with slope of 2 (see Fig b) rather than 
in exponential way in Levy walk model. 

Based on the above two facts, it is not a surprise to 
find that there are large deviations of the numerical 8 
values from the formula 8 = a — 1 (dotted line), see Fig 

E3 

In fact, the large deviation from the model used in 
Ref.|55j may result from inapproriate use of Levy walk. 
Levy walk is characterized by two assumptive distribu- 
tion functions which by specifying them appropriately, 
one can find a broad application spectrum for this sta- 
tistical model. However, as is stressed by Zaslavsky |3l) . 
great care must be taken in specifying these assumptive 
distribution functions (such as ipf(x, t) mentioned above) 
so that they do capture the characteristics of the dynam- 
ics. Failing to do so might be one of the reasons why 
8 = a — 1 obtained in Ref. |55j| can not be applied to 
the gas channels. Another reason is, in Levy walk it is 
assumed that between two successive steps there is no 
correlations; this is in clear contrast with the gas chan- 
nels. This leads us to believe that by using the Levy walk 
in an appropriate way (which the main dynamic features 
are incorporated and the correlations are well dealt with) 
one can get deeper insight of the transport behaviors in 
the gas channels and lattices. 



VII. CONCLUSIONS 

In summary, we have established a connection between 
anomalous heat conduction and anomalous diffusion in 
ID systems. The two exponents in anomalous difus- 
sion and anomalous thermal conductivity is connected 
via formula ©. This formula, includes all possible cases 
observed in different classes of ID models, ranging from 
subdiffusion, normal diffusion, and superdiffusion to bal- 
listic transport. Several conclusions can be drawn: (1) 
A normal diffusion leads to a normal heat conduction 
obeying the Fourier law. (2) A ballistic transport leads 
to an anomalous heat conduction with a divergent ther- 
mal conductivity kocL. (3) A superdiffusion leads to an 
anomalous heat conduction with a divergent thermal con- 
ductivity in thermodynamic limit. (4) More importantly, 
our result predicts that a subdiffusion system will be a 
thermal insulator. Existing numerical data from different 
low dimensional systems such as one dimensional lattice, 
small radius carbon nanotube and billiard gas channels 
support our results. 
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